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Summary 

The use of ceramics components in gas turbines, cutting 
tools, and heat exchangers has been limited by the relatively 
low flaw tolerance of monolithic ceramics. The development 
of whisker-toughened ceramic composites offers the potential 
for considerable improvement in fracture toughness as well 
as strength. However, the variability of strength is still too 
high for the application of deterministic design approaches. 
This report reviews several phenomenological reliability 
theories proposed for this material system, and reports on the 
development of a public domain computer algorithm. This 
algorithm, when coupled with a general-purpose finite element 
program, predicts the fast fracture reliability of a structural 
component under multiaxial loading conditions. 

Introduction 

The potential advantages of ceramic matrix composites 
include increased fracture toughness, and creep and corrosion 
resistance at very high service temperatures. The primary 
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applications under consideration are advanced turbine engine 
components, cutting tool bits, heat exchangers, and aerospace 
components (specifically those of the National Aerospace 
Plane). Considering that these composites will be produced 
from nonstrategic materials, it is not surprising that concerted 
research efforts are under way both in the field of materials 
science (to advance processing techniques) and in the field of 
engineering mechanics (to develop design methodologies for 
these material systems). 

The material system of interest in this report is the whisker- 
toughened ceramic matrix composite. Analysis of components 
fabricated from this material requires a departure from the 
design philosophy (i.e., the factor of safety approach) prevalent 
in designing metallic structural components, which are more 
tolerant of flaws. Since failure of components fabricated from 
this material is governed by the scatter in strength, statistical 
design approaches must be used. The primary objective of this 
report is to review several phenomenological failure models and 
to report on the development of a public domain computer 
algorithm which, when coupled with a general-purpose finite 
element program, -predicts the fast fracture reliability of a 
structural component under multiaxial loading conditions. The 
present version of this algorithm has been given the acronym 
TCARES (Toughened Ceramics Analysis and Reliability 
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Evaluation of Structures) and is a direct offspring of the CARES 
(a.k.a. SCARE) program (ref. 1), which has found widespread 
use in the design of monolithic ceramic components. 

In addition to capturing the inherent scatter in strength, the 
reliability analysis of components fabricated from whisker- 
toughened ceramics must account for material symmetry 
imposed by whisker orientation. A noninteractive macroscopic 
model has been presented that accounts for the transversely 
isotropic material symmetry (S.F. Duffy and S.M. Arnold, 
Noninteractive Macroscopic Statistical Failure Theory for 
Whisker Reinforced Ceramic Composites, to be published in 
J. Compos. Mater.) often encountered in hot-pressed and 
injection-molded whisker-toughened ceramics. A similar 
model (ref. 2) has been proposed for whisker-toughened 
ceramics with orthotropic material symmetry. This continuum 
approach excludes any consideration of the microstructural 
events that involve interactions between individual whiskers 
and the matrix. Other authors have addressed fracture of 
ceramic matrix composites on a more local scale. A model 
based on probabilistic principles has been developed to 
compute an increased energy absorption during fracture due 
to whisker pull-out (ref. 3). The processes of crack deflection 
(ref. 4) and crack pinning (ref. 5) have also been addressed. 
The latter two approaches are founded in deterministic fracture 
mechanics. Since these crack mitigation processes strongly 
interact, it is difficult to experimentally detect or analytically 
predict the sequence of mechanisms leading to failure. 

A more feasible approach is to compute reliability in terms 
of macrovariables by using a continuum-based criterion. This 
underscores a fundamental difference that exists between the 
materials scientist and the engineer. The materials scientist 
focuses on mechanisms of failure at the microstructural level, 
and the engineer focuses on this issue at the component level. 
The failure models currently incorporated into the computer 
algorithm TCARES adopt the engineer’s viewpoint. This point 
of view implies that the material element under consideration 
is small enough to be homogeneous in stress and temperature, 
yet large enough to contain a sufficient number of whiskers 
such that the element is a statistically homogeneous continuum. 
This does not imply that the microscopic and macroscopic 
levels of focus are mutually exclusive. Indeed, a close rela- 
tionship must exist between the materials scientist and the 
design engineer so as to develop better failure models to facil- 
itate the use of ceramic materials in structural components. 


Noninteractive Reliability Models 


the chain. Defining / as the probability of failure of an 
individual link gives 


/= AV (1) 

where AV denotes an increment in volume and ^ is a failure 
function per unit volume of material. By taking r as the 
reliability of an individual link, then 

r = l-t AV (2) 

If the failure of an individual link is considered a statistical 
event, and if these events are assumed to be independent, then 
the reliability of the continuum, denoted as /?, is given as 



where N denotes the number of links and \p(: t, ) is the failure 
function per unit volume at position x t within the continuum. 
Lowercase Roman letter subscripts here and in the following 
expressions denote tensor indices with an implied range from 
1 to 3. Greek letter subscripts and uppercase italic letters are 
associated with products or summations with ranges that are 
explicit in each expression. Alternatively, the reliability of the 
continuum is given by the following expression 



where the integral within the bracket is referred to as the risk 
of rupture. 

Depending on fabrication, a whisker-toughened composite 
may have isotropic, transversely isotropic, or orthotropic 
material symmetry. The principle of independent action (PIA) 
would be an appropriate first approximation macroscopic 
theory for isotropic whisker composites. In this instance the 
failure function \j/ would depend only on stress or the principal 
invariants of stress; that is, 

t = ^{o X ,0 2 ,a 3 ) (5) 

where is the Cauchy stress tensor and <J\, o 2 , and a 3 are 
the associated principal stresses. 

However, for transversely isotropic whisker composites, the 
failure function must also reflect material symmetry. This 
requires that 


Here, a continuum is considered to be a chain composed 
of links connected in series. Therefore, the overall strength 
of the continuum is governed by the strength of its weakest 
link. It is further assumed that the events leading to failure 
of an individual link are not influenced by any other link in 


\P = \l/(Oij,ai) ( 6 ) 

where a, is a unit vector that identifies a local material 
orientation. This orientation, depicted in figure 1, is defined 
as the normal to the plane of isotropy. The sense of a { is 
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SPECIMEN SPECIMEN 

Figure 1.— Examples of transversely isotropic whisker-toughened ceramics. 

immaterial, and thus its influence is taken through the product 
aflj\ that is, 

\p = (7) 

Note that apj is a symmetric second-order tensor, the trace 
of which satisfies the identity = 1. Furthermore, the 
stress and local preferred direction may vary from point to 
point in the continuum. Thus equation (7) implies that the stress 
field and the unit vector field (i.e., cr^**) and dj(x k )) must 
be specified to define 

For orthotropic composites the failure function must also 
reflect the appropriate material symmetry. This requires that 

^ = '{ / (0 i j>aja j ,b i bj) ( 8 ) 

where a x (a different vector than the one used for transverse 
isotropy) and b x are unit vectors that identify local material 
orientations. These vectors are assumed to be orthogonal such 
that afti = 0. 

Since \p is a scalar valued function dependent on second- 
order tensors, the form of \j/ must remain invariant under 
proper orthogonal transformations. This requires the function 
to be insensitive to the global coordinate system used to define 
the stress tensor and material directions. Through the use of 
invariant theory, a finite set of invariants known as an integrity 
basis can be developed for the isotropic, transversely isotropic, 
and orthotropic material symmetries (table I). The individual 
invariants of each integrity basis can be likened to a basis 
vector that helps to span a particular vector space (e.g., the 
set of unit vectors that span the Cartesian space). A slightly 
different set of invariants that correspond to physical mech- 
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TABLE I.— NONZERO MEMBERS OF 
INTEGRITY BASIS CORRESPONDING 
TO MATERIAL SYMMETRY 


Material symmetry 

Invariant 

Isotropy 

/, = a u 


h = a ij°ji 


h - a ij°jk a ki 

Transverse isotropy 

/, = a,, 


h = 


h - °ij a jk°ki 


h = afljOji 


h - a i a j°jk°ki 

Orthotropy 

l \ = °ii 


h- °ij a ji 


h - a ij°jk a ki 


h = 


h ~ a i a j°jk°ki 


A> = bjbjOjt 


If = bjbj<jj k o ki 


anisms related to failure is constructed from each integrity 
basis. (See table II for a brief description of each invariant 
and fig. 2 for a graphical interpretation.) These invariants can 
be identified with a principal stress or a component of the stress 
tractions coincident with a material direction. 

The invariants used to form \p are assumed to act inde- 
pendently in producing failure such that has the following 
general form: 



where N = 3 and M = 1 for isotropy, N = 4 and M = 3 for 
transverse isotropy, and N = M = 5 for orthotropy (table HI). 
In association with each invariant, the a’s correspond to 
Weibull shape parameters and the jS’s correspond to Weibull 
scale parameters. A variety of test methods could be used to 
determine these parameters. One approach is to obtain the data 
associated with the normal stress tractions from fast fracture 
of simple bend test specimens, often referred to as modulus 
of rupture (MOR) bars. The Weibull parameters associated 
with shear tractions would be obtained from appropriate shear 
strength tests. It is further assumed that compressive principal 
stresses and compressive stresses associated with a material 
orientation do not contribute to failure. 


TCARES Algorithm 

The basic data requirements of TCARES (fig. 3) closely 
follow the structure of its parent code CARES. The algorithm 
requires the stress analysis from a general-purpose finite 
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TABLE II. -INVARIANTS ASSOCIATED WITH PHYSICAL MECHANISMS 
DIRECTLY RELATED TO FAILURE 


Material symmetry 

Invariants used in \J/ function 

Comments 

Isotropy 

ft = °t 

Principal stresses; functionally dependent on the first three 


(l = o 2 
f\ - a 3 

invariants of stress 

Transverse isotropy 

i=h 

- Normal stress component of stress traction associated with a t 


s. 

rs 

i 

w-, 

ii 

Shear stress component of stress traction associated with a,- 


/ 3 = V6(/, - U) + *° 

Maximum normal stress in plane of isotropy 


u = '^(/, - U) - R a 

Minimum normal stress in plane of isotropy 

Orthotropy 

ft = h 

Normal stress component of stress traction associated with a t 


h = Us ~ W 2 Y A 

Shear stress component of stress traction associated with a i 


h = k 

Normal stress component of stress traction associated with b { 


h = lh - (4) 2 1* 

Shear stress component of stress traction associated with 6,- 


1 

rf 

1 

*>c 

II 

Normal stress in direction defined by cross product of a { and 6, 


°r = h ' a )/ 2 - / 5 + <*)< / 4 ) 2 - (^)(/|) 2 + mw * 




(a) Transverse isotropy. 

(b) Orthotropy. 

Figure 2.— Invariants associated with transverse isotropy and orthotropy. 


element code. Currently, the preliminary version of TCARES 
is compatible with MSC/NASTRAN, although it is anticipated 
that future versions compatible with the MARC, ADINA, and 
ANSYS finite element codes will be available. The algorithm 
allows the user to specify temperature-dependent statistical 
material parameters for each material symmetry. Alternatively, 
the program has the capability to estimate statistical parameters 
from fracture data obtained from uniaxial tensile or flexural 
specimens. (Details of this capability can be found in ref. 6.) 
The preceding section on noninteractive theoretical models 
implies a volume flaw analytical approach. It is quite possible 
that the surface and volume of a structural component will fail 
because of distinctly different flaw populations. Accordingly 
the TCARES program has the capability to separately conduct 
surface and volume reliability analyses. The program produces 
as bulk output a summary of input from the finite element code, 
element statistical properties, element survival probabilities, 
and an overall component survival probability. 

TCARES requires certain information from the finite ele- 
ment structural analysis. This includes element volumes, nodal 
temperatures, centroidal or nodal stresses, element principal 
stresses (for PIA analysis), and element identification numbers. 
The current version of TCARES assumes that the nodal 
stresses from the finite element code are provided relative to 
the local material orientation for transversely isotropic and 
orthotropic materials. This precludes having to input material 
orientation vectors for each finite element. 

The TCARES user input requirements are grouped into three 
categories. (See table IV.) The first category, entitled master 
control input, defines control indices for stress and graphics 
output, the number of component materials, and information 
regarding the finite element code and mesh. The second 
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TABLE III.— FUNCTIONAL FORMS OF 
CORRESPONDING TO MATERIAL SYMMETRY 




Figure 3.— Data requirements of TCARES. 


category, entitled material control input, allows the user to 
designate the statistical parameters, and defines control indices 
for analysis (volume or surface) and material symmetry (iso- 
tropic, transversely isotropic, or orthotropic). The final 
category contains the temperature-dependent material sym- 
metry statistical parameters. This includes the Weibull shape 
and scale parameters at each designated temperature. 

The output from a reliability analysis conducted with 
TCARES can be grouped into eight broad categories. They 
include an echo of the master control input, a summary of 
element types, stress output, volume, and averaged temperature 
of each element. An echo of the material control indices and 
a summary of the analysis of each element are also provided. 
The latter includes material identification, failure probability, 
risk of rupture intensity, and the Weibull statistical parameters 
for each element. The program sorts and identifies 15 elements 
with the highest risk of rupture intensity to aid the design 
engineer in locating “hot” areas within a component. Finally, 
the overall component survival probability is given. 

From Simple to Complex Geometries 

Multiaxial experiments are a necessity to assess the accuracy 
of the noninteractive modeling approach. One experimental 
avenue highlighted here is the application of thermomechanical 
loads to tubular specimens. Initially, a thick-walled tube 
subjected to an applied torque is considered. A second problem 
is presented where the same thick-walled tube is subjected to 
a simultaneous application of internal pressure and axial 
torque. In all applications considered, isothermal conditions 
are assumed. However, the algorithm is capable of non- 
isothermal analyses if the user specifies the values of the 
Weibull parameters at a sufficient (and appropriate) number 
of temperature values. Unfortunately, at the present time no 
data base exists to properly characterize the multiaxial 
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TABLE IV. -DESCRIPTION OF INPUT PARAMETERS 


Category 

Input 

Description 

Master control input 

NE 

Analysis option 

0— Experimental data analysis only 
1 — MSC/NASTRAN analysis 

2— ANSYS analysis (option currently unavailable) 

3— MARC analysis (option currently unavailable) 


NMATS 

Number of materials for surface flaw analysis 


NMATV 

Number of materials for volume flaw analysis 


IPRNT 

Control index for stress and/or fracture data output 

0— Suppress element stress and/or fracture data output 

1— Print element stress and/or fracture data output 


IGRPH 

Control index for graphics output 

0— Suppress element risk of rupture intensity output 

1— Print element risk of rupture intensity output 


NS 

Number of cyclic symmetry segments 

Temperature-independent 

MATID 

User-designated material identification number 

material control input 

ID1 

Control index for experimental data analysis 

1— Uniaxial tensile specimen data 

2— Four point bend test data 

3— All Weibull parameters specified by user 

4— Censored data for suspended item analysis 

(uniaxial tensile specimen) 

5— Censored data for suspended item analysis 

(four point bend test) 


ID2 

Control index for material symmetry failure criteria 

1— Isotropic (PI A) 

2— Transversely isotropic 

3— Orthotropic 


ID4 

ControLindex for- volume or surface flaw analysis 

1— Volume 

2— Surface 

T emperature-dependent 

TDEG 

Temperature for specified data set 

material control input 

PARAM 

Weibull parameters corresponding to TDEG 


statistical parameters for a whisker-toughened ceramic, 
although efforts (ref. 7) are under way to accomplish this goal. 
Thus, an assessment of the program output relative to actual 
structural component data is reserved for a later date. For the 
examples that follow, Weibull statistical parameters are 
assumed for the purpose of illustration; however, the values 
adopted are well within the range of the sparse data that can 
be found in the open literature (refs. 8 to 10). 

To test the validity of the reliability calculations performed 
by the program, a comparison of output with a hand calculation 
is presented for the aforementioned simple structural problem, 
that is, a thick- walled tube subjected to an axial moment or 
torque. It is assumed that the cylinder is fabricated from a 
whisker-toughened ceramic material having an orthotropic 
material symmetry such that a { = (0,0,1) and b { = (0,1,0) at 
every point in the structure. A cylindrical coordinate system 
readily lends itself to this application; hence, a, is directed 
along the z-axis of the cylinder and b t is oriented in the 0 


(circumferential) direction. With this geometry, material 
symmetry, and load condition, only the two terms (see 
table III) in the failure function associated with shear tractions 
are nonzero, and \ p takes the form 



where T is the applied torque, r is the radius, and J is the polar 
moment of inertia. Assuming an inner radius of 1 cm, an outer 
radius of 5 cm, and a length of 5 cm gives 


/? = exp^-jjj i pr dr dd dz 




(ID 
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Figure 4.— Risk of rupture intensity for thick-walled tube subjected to torque 
of 7500 N m. 

For dimensionless reliability, the Weibull scale parameters 
/J 5 ) have units of stress x (volume) 1/a , and the 
Weibull shape parameters are unitless. With a 2 = 10, 
0 2 = 15 000, a 4 = 6.5, j3 4 = 10 000 (the other Weibull 
parameters can be stipulated arbitrarily), and T = 7500 N m, 
the above integration yields an overall reliability of 83.6 
percent. The problem was also modeled by using MSC/ 
NASTRAN to generate a numerical solution for the stress 
distribution. The model was composed of 1500 eight-node 
elements (the stresses were within 2 percent of the closed form 
solution) which generated an overall component reliability of 
81.8 percent. Figure 4 is a color plot of the variation of \ p 
in a quarter section of the component. Note that ^ (risk of 
rupture intensity) is a measure of reliability independent of 
the element geometry, and that it attains a maximum value 
along the outer edge of the tube. 

Next, consider the same tube, subject to the conditions stated 
in the preceding paragraph, with an additional applied internal 
pressure of 70 MPa. Here a 3 and jS 3 can no longer be 
stipulated arbitrarily and take values of a 3 = 7.5 and 
/? 3 = 12 000. The overall component reliability decreases to 
77.4 percent. For this load case the circumferential stress is 
a maximum at the inner radius, and it decreases nonlinearly 
through the thickness. The shear stress from the applied torque 
is a minimum at the inner radius, and it increases linearly with 
the radius. Given this multiaxial stress distribution, one expects 
the maximum risk of rupture intensity to occur at some point 
midway through the thickness. However, this is not the case, 
as is evident in figure 5. The maximum risk of rupture intensity 


Figure 5.— Risk of rupture intensity for thick- walled tube subjected to torque 
of 7500 N m and internal pressure of 70 MPa. 

RISK OF 
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Figure 6.— Risk of rupture intensity for thick-walled tube subjected to torque 
of 7500 N m and internal pressure of 100 MPa. 


occurs at the inner radius, and much of the inner volume of 
the tube remains relatively “cold.” This underscores the need 
of not only considering overall component reliability, but also 
giving consideration as to where local “hot” spots occur within 
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a component. If this particular component were to fail, one 
would expect the failure to originate in the vicinity of the inner 
radius. 

Overall component reliability can be adversely affected by 
either increasing the stress distribution in additional regions 
of the component (the so-called size effect), or by dramatically 
increasing the stress locally (thereby increasing the chance of 
failing a single link in the chain). Increasing the internal 
pressure of the previous example to 100 MPa sharply decreases 
the component reliability to 36 percent. Figure 6 depicts the 
variation of \J/ throughout the component. It appears that the 
additional affected region of the component is minimal. How- 
ever, \p changes two orders of magnitude along the inner 
radius. Although not depicted here, the overall component 
reliability of the tube subjected to only an internal pressure 
of 100 Mpa is 44 percent, indicating that the major source 
in the degradation of reliability is the internal pressure. 

Concluding Remarks 

The applications presented here represent very straight- 
forward structural analyses. Similar calculations and graphical 
interpretations can be carried over into designs with much 
more complex geometry and boundary conditions. 

Recent advances in processing whisker-toughened ceramics 
have resulted in the reduction of inhomogeneities, uniform 
whisker distributions, and increasingly dense matrices, all of 
which have greatly improved the reliability of this material 
system. However, the variability of strength is still too high 
for the application of deterministic design approaches. Sta- 
tistical design methodologies must be used not only to account 
for the scatter in ultimate strength, but also to account for 
decreasing bulk strength with increasing component volume 
(the so-called size effect). If the orientation of the whiskers 
is such that an anisotropic material symmetry is imparted, this 


must also be accounted for. Several phenomenological failure 
theories that take into consideration these issues in a macro- 
scopic sense have been reviewed. In addition, a computer 
algorithm has been discussed that incorporates these theories 
and that is capable of predicting reliability given the state of 
stress and temperature distribution within a component. 
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